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摘 要 提出 一 种 基于 坐标 系 转换 与 QTM 编码 的 天 球面 探测 履 善 栅 格 化 分 析 方法 , 并 对 其 计算 精度 和 效率 进行 
了 实验 分 析 , 相 比 传统 算法 和 经 纬 网 格 编码 方法 , 该 方法 具有 明显 的 效率 和 精度 优势 , 能 够 满足 以 天 球面 为 探测 目 
标的 空间 科学 卫星 有 效 载荷 探测 覆盖 分 析 和 任务 规划 的 需要 
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Abstract In sky coverage analysis, we present an efficient algorithm based on a coordinate trans- 
formation method to calculate the points on the boundary of detection field. In this method, the 
complexity of coordinates transformation is reduced by building a new Celestial Coordinates. Based 
on this algorithm, QTM, a coding method that has many advantages can be used to manage the 
data calculated. The main contents of this topic include: to build the quantitative calculation model 
for sky coverage; to establish the rasterized expression method of sky coverage analysis results; and 
to verify the method, analyse the accuracy of the method, compare the efficiency with the attitude- 
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0 引言 


中 国 科学 院 空间 科学 先导 专项 计划 在 “十 二 五 ” 
期 间 发 射 5 颗 空 间 科学 卫星 , 其 中 , 硬 X 射线 调制 户 
远 镜 卫 星 和 暗物质 粒子 探测 卫星 的 探测 对 象 是 硬 X 
射线 、 高 能 电子 和 伽 马 射线 等 宇宙 辐射 , 将 采用 巡天 
观测 、 定 点 观测 等 不 同 观测 模式 对 空间 进行 观测 . 在 
对 这 两 颗 卫 星 进行 地 面 运行 控制 的 过 程 中 , 需要 频繁 
计算 一 定时 间 段 内 其 探测 载荷 在 不 同位 置 、 姿 态 和 
工作 模式 下 的 探测 视 场 对 探测 目标 或 探测 区 域 的 覆 
盖 情 况 , 作为 任务 规划 的 输入 条 件 和 依据 , 这 是 通过 
任务 规划 方法 寻求 最 优 控制 方案 的 前 提 条 件 趾 . 

天 文 观 测 等 对 天 观测 的 观测 对 象 和 探测 覆盖 范 
围 一 般 采 用 天 球 坐标 系 来 描述 , 天 球 坐 标 系 是 天 文学 
中 常用 的 坐标 系 四 ， 在 天 文学 中 , 天 球 是 一 个 假想 
的 球 , 理论 半径 无 限 大 , 球 心 即 地 心 . 天 空中 所 有 物 
体 均 被 投影 到 天 球 的 球面 上 , 基于 地 心 赤道 惯性 坐标 
系 回 的 天 球 与 地 球 相对 应 , 也 具有 天 赤道 和 天 极 . 其 
中 天 赤道 即 地 赤道 在 天 球 上 的 投影 , 天 极 即 地 球 极点 
在 天 球 上 的 投影 . 天 球面 投影 如 图 1 所 示 . 星空 界面 
中 RA (Right Ascension) 为 赤 经 , 表示 天 体 投影 点 到 


本 


north 
celestial 


/equinox 
point 


图 1 天 球面 投影 示意 
Fig.1 Celestial projection 


Celestial surface, Cavalcade approach method, Coordination translation, QT™, 


春分 点 的 角度 距离 , 用 时 、 分 、 秒 来 表示 ; 偏差 (Dec- 
lination) 为 赤 纬 , 表示 天 体 投影 点 到 天 赤道 的 角度 距 
离 , 单位 用 度 、 分 、 秒 来 表示 四 

在 近 地 卫 星 的 天 球 而 探测 覆盖 分 析 中 , 对 于 星体 
而 言 , 卫星 位 置 的 变化 可 以 忽略 不 计 , 其 探测 视 场 在 
天 球面 上 投影 的 形状 和 大 小 保持 不 变 , 因而 只 需求 解 
出 观测 矢量 在 天 球面 上 投影 点 的 赤 经 和 未 纬 即 可 求 
解 天 球面 瞬时 覆盖 范围 

但 是 , 当前 中 国 在 天 球面 覆盖 分 析 方面 的 关注 和 
研究 较 少 , 现 有 成 熟 的 对 地 探测 覆盖 计算 模型 在 应 用 
于 天 球面 覆盖 分 析 时 存在 一 些 局 限 性 , 在 效率 及 计算 
精度 等 方面 难以 满足 天 球面 探测 覆盖 快速 计算 的 要 
求 , 需要 开展 专门 的 柳 盖 分 析 方 法 研究 . 设计 实现 天 
球面 探测 覆盖 分 析 方法 , 需 重点 解决 的 技术 问题 主要 
包括 以 下 两 个 方面 

(D) 计算 效率 

在 卫星 运行 控制 过 程 中 , 一 次 任务 规划 涉及 的 时 
间 周期 非常 长 (一 个 月 甚至 一 年 ) 在 如 此 大 的 时 间 路 
度 内 , 需 按照 一 定时 间 步 长 根据 地 球 公转 、 卫星 轨 首 
位 置 、 姿 态 参 数 和 工作 模式 等 进行 频繁 计算 , 运算 量 
巨大 ， 如 何在 确保 计算 精度 的 前 提 下 通过 算法 改进 
提高 计算 效率 和 速度 , 是 工程 实现 中 需要 重点 解决 的 
问题 

(2) 覆盖 计算 结果 的 机 格 化 表达 

由 于 有 效 载荷 探测 视 场 一 般 采 用 解析 式 的 几何 
形状 来 描述 , 而 任务 规划 算法 需要 以 棚 格 化 的 覆盖 计 
算 结果 作为 输入 , 为 实现 对 天 球面 剖 分 编码 效率 、 数 
据 处 理 精 度 及 计算 速度 等 各 个 指标 的 兼顾 , 应 考虑 选 
用 合适 的 球面 网 格 化 方法 来 管理 探测 目标 和 探测 区 
域 , 并 实现 计算 结果 的 快速 转换 


1 基于 坐标 系 转换 的 天 球面 探测 
履 资 快速 算法 
1.1 ”空间 载荷 探测 视 场 的 数学 描述 


对 空间 进行 观测 的 卫星 有 效 载荷 探测 视 场 , 按 其 
空间 属性 可 以 划分 为 以 下 两 类 回 . 
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(1) 矩形 猴 视 场 

其 视 场 形状 可 使 用 两 个 方向 的 张 角 a, 2 来 描述 ， 
矩形 锥 探测 视 场 的 4 条 边 在 天 球面 上 的 投影 是 4 个 
顶点 4, B,C,D 之 间 的 球面 大 圆 弧 , 有 效 载荷 矩形 锥 
视 场 在 天 球面 上 的 探测 覆盖 范围 是 由 这 4 条 球面 大 
圆 弧 所 围 成 的 区 域 , 其 中 4B, CD 球面 大 圆 劣 弧 长 
为 5, 而 4D, BC 球面 大 圆 劣 弧 长 为 % 如 图 2 所 示 . 

(2) 圆锥 视 场 

其 视 场 形状 可 使 用 圆锥 的 轴 及 半 张 角 入 来 描述 ， 
圆锥 探测 视 场 的 轴 与 天 球面 的 交点 为 P, 圆锥 探测 
视 场 的 各 条 母线 与 天 球面 的 交点 形成 一 个 以 已 为 中 
心 的 球面 小 圆 , 有 效 载 荷 矩形 锥 视 场 在 天 球面 上 的 控 
测 覆 盖 范 围 是 由 以 天 球面 已 点 为 圆心 、 半 径 为 ^ 的 
球面 小 圆 所 围 成 的 区 域 , 如 图 3 所 示 . 
1.2 ”快速 求解 算法 

由 空间 载荷 探测 视 场 的 数学 描述 分 析 可 知 , 其 覆 
盖 求 解 可 简化 为 球面 大 圆 弧 段 的 求解 问题 和 已 知 圆 
心 与 球面 半径 的 球面 小 圆 求解 问题 ， 关 于 球面 大 圆 
的 求解 问题 , 在 航空 、 航 海 等 领域 已 有 一 些 成 卉 算法 . 

洪 德 本 提出 了 一 种 用 解析 法 设计 大 圆 航线 的 计 
算 方法 , 并 将 其 用 于 计算 机 解 算 名, 通过 解 算 大 圆 航 
线 , 求 出 大 圆 航路 上 各 点 坐标 、 航 向 及 航程 . 该 算法 
利用 球面 三 角 余 弦 定 理 计 算得 到 大 圆 航程 , 然后 将 大 
圆 航程 按 一 定 间隔 划分 为 若干 段 , 各 段 的 距离 即 为 各 
航路 点 之 间 的 航程 , 由 球面 三 角 正 弱 定 理 可 得 各 段 分 
割 点 ( 即 航路 点 ) 的 坐标 . 该 算法 的 计算 比较 简单 , 并 
且 可 以 根据 计算 精度 需要 调整 各 分 割 点 之 间 的 间隔 . 
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图 2 甜 形 锥 视 场 探测 覆盖 范围 
Fig.2 Projection to the celestial surface of rectangle 


awl detection field 
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但 是 航线 设计 的 主要 目的 是 确定 航向 角 , 因而 每 个 分 
割 点 均 须 计算 航向 角 , 增加 了 计算 量 , 并 且 由 于 其 近 
似 认为 在 每 个 分 割 弧 段 内 的 航向 角 为 恒定 值 , 造成 了 
累积 误差 , 不 能 够 适应 天 球面 探测 覆盖 计算 的 需要 . 
对 于 球面 小 圆 的 求解 , 通常 采用 在 三 维 直 角 坐 标 系 下 
求解 出 球面 小 圆 各 点 的 坐标 值 , 再 将 其 转换 为 经 纬度 
坐标 的 方法 . 

针对 上 述 算法 的 不 足 , 本 文 提出 一 种 基于 坐标 系 
转换 的 天 球面 探测 覆盖 快速 算法 , 在 天 球 坐 标 系 下 取 
任意 一 点 N' (其 在 天 球 坐 标 系 下 的 坐标 为 赤 纬 和 0， 
赤 经 po) 作为 新 的 天 极点 , 以 通过 N' 的 天 球面 大 圆 
弧 为 0"” 赤 经 , 建立 新 的 球面 坐标 系 , 新 建立 的 球面 
坐标 系 与 天 球 坐 标 系 的 几何 关系 如 图 4 所 示 . 
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图 3 圆锥 视 场 探测 覆盖 范围 


Fig.3 Projection to the celestial surface of circular 


cone detection field 


图 4 新 建 球面 坐标 系 与 原 天 球 坐标 系 的 几何 关系 
Fig.4 ， Relationship between the two celestial 


coordinates 
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新 建立 的 球面 坐标 系 与 天 球 坐标 系 类 似 , 也 是 
球面 坐标 系 , 不 同 的 是 新 建 球面 坐标 系 的 2 轴 指 向 
是 N 点 方向 , (Xo,po) 是 N' 点 的 天 球 坐 标 , 与 其 垂 
直 的 大 圆 面 为 球面 坐标 系 的 赤道 , X' 轴 和 Y' 轴 在 
其 赤道 平面 内 , N' 点 和 北 天 极点 N 确定 的 子午 线 为 
零度 子午 线 . 天 球面 的 一 点 由 纬度 X (该 点 与 球 心 0 
点 连 线 与 赤道 的 夹 角 ) 和 经 度 (过 该 点 的 子午 面 与 
过 北 天 极点 的 子午 面 之 间 的 夹 角 ) 来 标识 . 

对 于 天 球面 上 任意 一 点 , 其 在 新 建立 的 球面 坐标 
系 下 的 坐标 (入, yp") 与 天 球 坐 标 系 下 的 坐标 (和 ,yp) 之 
间 的 转换 关系 为 门 : 


coS 入 = — cosAcosAo— 


sin 入 Sin Mo cos(% 一 20)， 


sinw' 一 Sin 和 sin(Pp 一 po)csc 和 /， (1) 
cos 入 一 cos Acos X0 一 Sin X sin Mo cos P/， 


sin(p 一 %0) = SinAX sin %'csc 入 . 


根据 探测 视 场 边界 的 几何 特点 , 建立 新 的 天 球面 
坐标 系 , 在 新 的 坐标 系 下 进行 求解 , 可 以 大 幅度 减少 
计算 量 . 
1.3 ”算法 实现 

对 于 和 矩形 锥 探测 视 场 , 选取 矩形 视 场 一 条 边 的 一 
个 顶点 为 新 球面 坐标 系 下 的 北极 点 , 定义 矩形 视 场 该 
条 边 所 在 的 半 贺 弧 为 新 球面 坐标 系 下 的 0” 赤 经 , 则 
矩形 锥 视 场 该 条 边 在 新 球面 坐标 系 下 经 度 为 0", 纬 
度 范围 在 (90° 一 a) 或 (90° 一 6) 至 90° 之 间 , 然后 依 
据 式 (1) 可 将 矩形 视 场 该 条 边 上 各 点 在 天 球面 坐标 
系 下 的 赤 经 未 纬 坐标 快速 求 出 , 如 图 5 所 示 . 

对 于 圆锥 探测 视 场 , 选取 视 场 中 心 点 为 新 球面 坐 
标 系 下 的 北极 点 , 则 圆锥 视 场 边界 线 在 新 球面 坐标 系 
下 的 经 度 在 0" 至 360° 之 间 , 纬度 为 (90° - 入, 然后 
同样 依据 式 (1) 可 将 圆锥 视 场 边界 上 各 点 在 天 球面 
坐标 系 下 的 赤 经 赤 纬 坐标 快速 求 出 , 如 图 6 所 示 . 在 
算法 的 实现 中 , 式 (1) 中 的 正弦 函数 在 第 一 、 二 象限 
均 为 正 , 有 时 会 出 现 两 解 . 可 以 使 用 球面 四 联 公 式 代 
替 正 弦 公 式 进 行 求解 运算 , 消除 解 的 不 确定 性 并 提高 
计算 精度 四. 
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图 5 甜 形 锥 视 场 探测 覆盖 范围 快速 算法 示意 
Fig.5 Sketch of the calculation of the rectangle 


awl detection field 
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图 6 圆锥 视 场 探测 覆盖 范围 快速 算法 示意 
Fig.6 Sketch of the calculation of the circular 


cone detection field 


图 7 正八 面体 的 构造 与 编码 
Fig.7 Establish and code of the regular octahedron 
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2 基于 QTM 编码 的 天 球面 
探测 覆盖 栅 格 化 表达 


球面 网 格 系统 是 一 种 等 积 投影 系统 ， 其 精度 可 
随 网 格 精细 程度 增加 而 增加 ， 适 合 于 探测 覆盖 情况 
的 分 析 和 探测 任务 规划 . 在 众多 球面 网 格 化 算法 中 ， 
Duttong-I 提出 的 基于 正八 面体 的 球面 四 元 三 角 
网 (Quaternary Triangular Mesh, QTM), 因 其 层次 
性 结构 ,成 为 目前 研究 全 球 多 尺度 海量 数据 结构 的 有 
效 管理 方法 之 一 [9. 
2.1 QTM 构造 与 编码 

QTM 基于 正八 面体 进行 构造 , 正八 面体 的 上 下 
两 顶点 与 南北 天 极点 重合 , 其 棱 分 别 与 0", 90°, 180?， 
270° 赤 经 , 0" 赤 纬 的 投影 重合 03. 如 图 7 所 示 , 其 
对 应 的 编码 为 0~7. 

随后 , 逐 层 对 三 角 面 进行 编码 , 每 一 次 细 分 中 , 一 
个 三 角 面 被 分 成 4 个 新 的 三 角 面 , 中 间 三 角形 为 0， 
左面 三 角形 为 2, 右面 三 角形 为 3 上 面 或 下 面 的 三 
角形 为 1, 如 图 8 所 示 . 

一 个 上 层 的 QTM 编码 由 aoaiazas .… ak 表示 ， 
其 中 ao 为 一 个 八 分 码 , 表示 最 初 的 正八 面体 的 某 一 
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个 面 , a14a24a3…ax 为 个 四 分 码 , 表示 第 层 的 三 
角 面 片 编码 ,上 = 6 时 的 QTM 网 络 划分 结果 如 图 9 
所 示 . 
2.2 ” 赤 经 赤 纬 到 QTM 编码 的 转换 算法 

本 文 基于 赵 学 胜 等 提出 的 行列 通 近 法 (Caval- 
cade Approach Method)(13 进行 赤 经 赤 纬 到 QTM 码 
的 快速 转换 . 这 是 一 种 球面 小 圆 的 剖 分 方法 . 其 基本 
方法 为 , 先 根 据 剖 分 层次 求 出 QTM 的 最 大 行列 数 ， 
再 根据 所 求 点 的 经 纬度 计算 对 应 的 三 角 格 网 所 在 行 
列 数 , 然后 逐 层 递 推 求 出 对 应 的 QTM 编码 . 
3 实验 结果 分 析 

基于 所 提出 的 坐标 系 转换 与 QTM 编码 方法 , 本 
文 开 发 了 一 套 天 球面 探测 覆盖 栅 格 化 分 析 软 件 , 对 分 
析 结 果 的 正确 性 进行 验证 , 如 图 10 所 示 , 并 对 该 方法 
的 精度 和 效率 进行 了 分 析 . 
3.1 ”计算 效率 分 析 

计算 中 所 使 用 的 实验 平台 为 台式 计算 机 . CPU 
型 号 为 mtel@ CoreIM 2 Duo E8400, 主 频 3.0 GHz; 
内 存 4G RAM; 操作 系统 为 Windows XP. 


图 8 细 分 过 程 及 其 编码 
Fig.8 Subdivisions and its code 


图 9 一 个 6 层 的 QTM 网 格 
Fig.9 A 6-layer QTM mesh 


表 1 计算 1000 个 时 相片 段 的 矩形 锥 视 场 的 效率 对 比 
Table 1 Comparison in efficiency of rectangle 


awl detection field 


分 辨 率 大圆 航 迹 姿态 矩阵 转换 。 球面 坐标 系 
/(") ”法 耗 时 /ms ”法 耗 时 /ms ”转换 法 耗 时 /ms 
36 2146 1166 827 
18 4289 2334 1629 
个 10750 5792 4098 
3.6 21550 11588 8208 
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(a) 


图 10 8 层 编码 下 天 球面 覆盖 计算 结果 示例 . (a) 圆锥 视 场 (b) 矩形 视 场 


Fig.10 Coverage in an 8-layer-QTM mesh based on the calculation: (a) circular cone detection field, (b) rectangle 


awl detection field 


表 2 计算 1000 个 时 相片 段 的 圆锥 视 场 的 效率 对 比 
Table 2 Comparison in efficiency of circular 


cone detection field 


分 辩 率 /(“) ”姿态 年 阵 转换 球面 坐标 系 
法 耗 时 /ms 转换 法 耗 时 /ms 
36 1371 845 
18 2773 1678 
人 6840 4163 
3.6 13610 9310 


对 于 矩形 锥 视 场 的 计算 效率 , 这 里 采用 球面 坐 
标 系 转换 算法 与 大 圆 航线 解析 法 进行 效率 对 比 . 计 
算 1000 个 时 相片 段 , 结果 列 于 表 1 (分 辨 率 基于 张 角 
为 10° 的 情况 ). 可 以 看 出 , 在 计算 1000 个 时 相片 段 
时 , 大 圆 航 迹 法 每 次 都 要 计算 航向 , 效率 最 低 . 基于 坐 
标 系 转换 法 时 , 须 先 计算 出 视 场 在 卫星 本 体 坐 标 系 下 
的 覆盖 情况 , 利用 覆盖 随 卫星 共同 运动 的 特点 , 通过 
姿态 矩阵 计算 天 球面 坐标 系 下 的 覆盖 情况 , 效率 大 幅 
提升 . 二 是 利用 球面 坐标 系 , 相 比 较 于 姿态 算 阵 , 是 
二 维 坐标 到 二 维 坐标 的 转换 , 其 效率 也 有 一 定 提升 . 

对 于 圆锥 视 场 , 通常 都 是 使 用 姿态 矩阵 转换 法 这 
一 相对 通用 的 计算 方法 .但 是 姿态 矩阵 法 没有 考虑 
视 场 本 身 所 具有 的 一 些 特殊 性 质 . 对 于 圆锥 视 场 , 其 
在 新 建 球面 坐标 系 下 为 一 赤 纬 圈 , 具有 完全 对 称 的 特 
性 , 在 计算 时 , 利用 对 称 性 , 仅 需 计 算 0” ~180° 的 范 
围 即 可 . 效率 对 比 结果 列 于 表 2 (分 辩 率 基于 圆锥 张 
角 10* 的 情况 ). 
3.2” QTM 编码 精度 

QTM 中 某 一 级 编码 对 应 的 三 角 格 网 边 长 可 代 
表 此 级 别 的 分 辨 率 [1 站. 在 天 球面 覆盖 分 析 中 , 可 用 


表 3 ”整体 计算 效率 
Table 3 Efficiency in all 
重复 计算 次 数 ”矩形 锥 耗 时 /s ”圆锥 耗 时 /s 
1000 1.071 0.867 
5000 5.364 4.372 
10 000 10.811 8.686 


三 角 格 网 边 长 对 应 的 球 心 角 表示 当前 编码 级 别 的 分 
辩 率 . 本 文 基于 快速 转换 算法 , 利用 经 纬度 平分 法 进 
行 球 面 三 角形 的 细 分 , 其 三 角形 在 大 小 和 形状 上 不 可 
避免 地 有 一 些 形变 . 以 0? 赤 纬 所 在 边 作为 分 状 率 标 
准 , 对 于 大 层 的 QTM 编码 , 其 在 0? 赤 纬 上 有 2* 条 
边 , 则 对 应 的 球 心 角 为 90°.2*, 即 对 于 大 层 的 QTM 
其 分 状 率 为 90" . 2 下 . 

在 行列 台 近 法 的 算法 中 , 因 其 只 计算 了 奇数 列 ， 
故 存在 额外 的 半 个 格 网 的 误差 03]， 综 合 分 辩 率 约 
为 1.5 x 90°? . 2-*， 当 QTM 的 编码 层 数 有 取 19 
时 , 其 综合 分 辨 率 可 达到 角 秒 量 级 . 
3.3” 栅 格 化 覆盖 分 析 的 综合 精度 与 效率 

在 实际 工程 应 用 中 , QTM 的 编码 层 数 选取 以 满 
足 履 盖 分 析 的 实际 精度 需要 为 准 , 例如 当 卫 星 视 场 张 
角 达 到 30° 时 , 取 QTM 的 编码 层 数 = 12, 对 应 的 
栅 格 化 覆盖 分 析 精 度 为 0.03°, 能 够 满足 日 常任 务 分 
析 和 规划 的 需求 , 其 整体 效率 列 于 表 3. 


4 结语 

实验 验证 结果 表明 , 本 文 所 提出 的 天 球面 探测 覆 
盖 栅 格 化 分 析 方 法 相 比 于 非 栅 格 化 的 履 盖 分 析 方 法 ， 
具有 计算 效率 高 , 精度 可 控 , 计算 结果 易于 存储 和 表 
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达 , 交 并 运算 算法 易于 实现 等 技术 优势 . 在 实际 工程 
应 用 中 , 将 分 布 在 天 球面 上 的 探测 目标 (区 域 目标 或 
点 目标 ) 按照 QTM 进行 编码 , 表达 为 一 个 或 一 组 编 
码 值 (一 个 球面 三 角形 或 者 一 系列 球面 三 角形 拼接 
得 到 的 区 域 ), 就 可 以 直接 通过 编码 检索 的 方法 从 探 
测 覆 盖 分 析 结 果 中 快速 查找 该 目标 的 探测 机 会 , 从 而 
很 好 地 满足 任务 规划 算法 的 输入 需求 . 


该 方法 在 实际 应 用 中 还 可 根据 工程 约束 条 件 , 在 


分 析 结 果 的 存储 管理 效率 等 方面 进行 进一步 优化 . 
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